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Мультиагентный метод оптимизации 
с адаптивными параметрами” 


Исследован мультиагентный подход к оптимизации многомерных нелинейных функций. Разработан 
мультиагентный метод оптимизации с адаптивными параметрами. Проведены эксперименты по поиску 
глобального оптимума многомерной функции на основе предложенного метода. 


Введение 


Градиентные методы безусловной оптимизации, традиционно применяемые при 
синтезе моделей сложных объектов и систем, являются высокоитеративными и накла- 
дывают определенные требования (например, унимодальность, непрерывность, моно- 
тонность, дифференцированность и т.п.) к виду целевой функции, а также, как правило, 
являются методами локального поиска. 

Поэтому актуальной является задача разработки новых методов оптимизации, не 
выдвигающих подобные требования к целевым функциям. В связи с этим возникает ин- 
терес к методам, основанным на случайном поиске, которые характеризуются большей 
скоростью работы, а также не накладывают никаких требований на зависимости, кото- 
рые исследуются. К таким методам относятся мультиагентные методы интеллектуальной 
оптимизации, которые имеют бионическую природу и основаны на моделировании по- 
ведения коллективного интеллекта социальных живых существ (5\уатт шеШеепсе) [1], [2]. 
К интеллектуальным методам мультиагентной оптимизации, основанным на моделирова- 
нии общественного интеллекта, относятся: метод муравьиных колоний (Ап Соопу Ор- 
Иплайоп, АСО) [3], [4]; метод пчелиной колонии (Вее Со]опу ОрНпитаноп, ВСО) [5], [6]; 
оптимизация с использованием роя частиц (Рагасе З\уапи ОрНпитайоп, РЗО) [7] и 
другие методы [8]. Данные методы уже эффективно применяются для решения раз- 
ных задач оптимизации [3-10]. 


Постановка задачи 


Метод оптимизации с использованием роя частиц традиционно используется для 
решения задачи поиска глобального оптимума многомерной функции, где показал до- 
вольно хорошие результаты [11], [12]. Однако РЗО-метод характеризуется и недостат- 
ками, связанными с возможностью слишком быстрого изменения скорости частиц, что 
приводит к увеличению времени поиска оптимального решения. 
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Цель данной работы -— разработка модификации метода оптимизации с исполь- 
зованием роя частиц путем введения дополнительных процедур с целью исключения его 
недостатков, связанных с лавиноподобным изменением скорости частиц, что позволит 
улучшить оптимизационный процесс. 

Основными задачами работы являются: 

— исследование метода оптимизации с использованием роя частиц; 

— выявление основных недостатков исследуемого метода; 

— создание метода оптимизации на основе моделирования поведения роя частиц, 
исключающего основные недостатки базового метода; 

— разработка программного обеспечения, реализующего предложенный метод, и 
его тестирование. 


Мультиагентный метод оптимизации 
с использованием роя частиц 


В Р5О-методе агенты, называемые частицами, перемещаются в многомерном про- 
странстве решений. Изменения координат частиц внутри пространства поиска обус- 
ловливаются естественной социально-психологической тенденцией частиц конкури- 
ровать между собой. Итак, изменения состояний частицы зависят от опыта и знаний ее 
соседей. В этом случае слово «знание» является синонимом «информации». Результат 
такого моделирования заключается в том, что процесс поиска мотивирует частицы 
недетерминированным образом возвращаться в оптимальные участки пространства 
решений. 

Агенты в рое частиц имеют очень простое поведение: они стремятся превзойти 
достижение соседних частиц и улучшить собственные. Таким образом, емерджентное 
свойство данной системы состоит в исследовании оптимальных участков многомер- 
ного пространства поиска. 

Р5О-метод управляет роем частиц и каждая частица представляет собой потен- 
циальное решение. По аналогии с эволюционными стратегиями, рой можно тракто- 
вать как популяцию, а частицу как индивида (или хромосому). Таким образом, частица 
перемещается в многомерном пространстве решений и ее позиция определяется исхо- 
дя из собственного опыта и опыта своих соседей. 

Метод РЗО может быть описан в виде последовательности таких этапов: 

— инициализация; 

— вычисление значений целевой функции; 

— обновление скоростей и позиций частиц; 

— проверка критериев останова поиска. 

Этап 1. Инициализация. 


На данном этапе создается множество агентов А (|А| = И, ), определяющих точ- 


ки в И, -мерном пространстве поиска. 


Через х{1) обозначим позицию частицы # в пространстве поиска в момент време- 
ни { (1 обозначает дискретные значения времени). Позиция частицы меняется путем 
добавления скорости у) к текущей позиции: 

хкЕ+ 1) =хК® + УЕ 1). 


Начальное состояние определяется таким образом: х,(0) =ИЦ(х ), где (С(а,Ь) 


является функцией генерации случайных чисел из диапазона [а, 6]. Данная формула 
представляет собой вектор скорости и определяет сам оптимизационный процесс, а 
также отображает использование как полученных знаний частицы, так и обмен инфор- 
мацией с соседними частицами. Собственные знания самой частицы, которые также 


шы ? Х тах 
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называются когнитивным компонентом формулы скорости, прямо пропорциональны 
текущему расстоянию частицы от ее наилучшего положения, которое было найдено 
с момента начала ее жизненного цикла. А обмен информацией данной особи (агента) 
с другими агентами является социальным компонентом формулы скорости. 

Этап 2. Вычисление значений целевой функции и определение наилучших позиций. 
Данный этап может быть выполнен с помощью следующих шагов. 

Шаг 1. Установить: 7 = 1, где {— счетчик, определяющий номер текущего агента. 


Шаг 2. Определить локально наилучшую позицию ),. Если выполнено условие 
Л(х,) < Г(»,) ‚ тогда установить: у; = Х,, где Х, — текущая позиция 1-й частицы; /(х,) 
и /(у,) — значения целевой функции {в точках Хх; и У, соответственно. 


Шаг 3. Определить глобально наилучшую позицию у*. Установить: у* = = 
Шаг 4. Установить: #=1+1. 
Шаг 5. Если 1 < И, ‚ тогда выполнить переход к шагу 1, в противном случае — пе- 


рейти к выполнению следующего этапа. 
Этап 3. Обновление скоростей и позиций частиц. 
Известны два подхода, которые являются разновидностями базового РЗО-метода: 
зЪея и ея, которые отличаются степенью связанности частиц в пространстве поиска. 
В 26е51 РЗО-методе скорость частицы рассчитывается по формуле: 


уу (+1 =у,(+ет (Оу - ху] сз КО - ху, 
где у,(1) — скорость частицы # в измерении / (/=1,..., п, ) в момент времени 1; х,(1) — 


позиция частицы # в измерении /; у*(1) — наилучшая достигнутая позиция роя (в про- 
странстве решений); с1 и с — положительные константы ускорения, которые исполь- 
зуются для варьирования весов когнитивной и социальной компонент скорости частицы 
соответственно; 7/,(1),7, (1) =( (0,1) являются случайными значениями из диапазона 
[0, 1]. Эти случайные величины вносят стохастический элемент в работу метода. 

Величина у; отображает наилучшую позицию частицы 1, которую она посещала, 
начиная с первой итерации. Следующая оптимальная позиция частицы 7 в момент вре- 
мени {+ 1 рассчитывается по формуле: 


УК), если Л(х +1) > Л(р(Ю)}; 
х, (1+ П,если Л(х,(1+1)) < Л(у:(0), 


где Г :53"х —> \ — целевая функция (фитнесс-функция), 33"х — множество значений не- 


УП) = 


зависимых переменных, 5} — множество значений оптимизируемой величины. Так же, 

как и в эволюционных подходах, фитнесс-функция является мерой близости данного 

решения к оптимальному, а также определяет производительность или качество частицы. 
В Без{ РЗО-методе скорость рассчитывается по формуле: 


у (1+9 =у, (чай (Оу @) — ху] + со (Оу *, @-х,@} 
где у*, (1) — наилучшая позиция, найденная соседями 1-й частицы в измерении /. Ло- 
кально наилучшая позиция частицы у*, , то есть лучшая позиция, найденная в соседст- 
ве М,, которое определяется как: 


у*, (+1 е М, | 1(*, @+1)} = ши Ух), Ухе М}, 
где соседство М, определяется по формуле: 


О 


для соседства размерностью их. 
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Алгоритм этапа обновления скоростей и позиций частиц РЗО-метода может быть 
представлен в виде следующей последовательности шагов. 

Шаг 1. Установить: 1 = 1. 

Шаг 2. Обновить скорость частицы, используя формулу 25е3{ РЗО или е$Е РЗО. 

Шаг 3. Обновить позицию частицы. 

Шаг 4. Установить: #=1+1. 

Шаг 5. Если 1< п, ‚ тогда выполнить переход к шагу 2, в противоположном слу- 


чае — перейти к этапу проверки критериев окончания поиска. 

Этап 4. Проверка критериев окончания поиска. 

Окончание выполнения Р5О-метода происходит, если выполняется хотя бы один 
из критериев останова. 

Существуют два аспекта, которые должны быть учтены при выборе подобных мер. 

1. Условие останова не должно служить причиной преждевременной сходимости 
Р5О-метода. В противоположном случае будут получены субоптимальные решения. 

2. При проверке критериев остановки фитнесс-функция должна вызываться ми- 
нимальное количество раз, поскольку от этого сильно зависит скорость работы Р$ЗО- 
метода. 

В качестве критериев останова, как правило, используют: 

— достижение максимально допустимого количества итераций; 

— нахождение приемлемого решения; 

— отсутствие существенного улучшения значения целевой функции на протяже- 
нии некоторого количества итераций; 

— выход при условии, что первая производная фитнесс-функции в окраине части- 
цы или частиц равняется нулю. 


Мультиагентный метод оптимизации 
с адаптивными параметрами 


Эффективный метод оптимизации должен обладать двумя важными свойствами: 
максимально полно исследовать пространство решений задачи и дать возможность де- 
тализации поиска возле потенциальных оптимумов. В РЗО-методе это определяется 
формулой скорости частиц. 

В ранних реализациях РЗО-метода было выявлено, что скорости могут слишком 
быстро возрастать, особенно скорости тех частиц, которые находятся далеко от своего 
собственного оптимума или общего оптимума своих соседей. В результате, такие час- 
тицы могут оставлять пространство поиска, что очень отрицательно влияет на работу 
метода в целом. 

Для предотвращения подобных ситуаций в разработанном мультиагентном ме- 
тоде оптимизации с адаптивными параметрами предложено контролировать диапазон 
изменения скоростей частиц. Если скорость частицы превышает максимально допус- 
тимую, то ее следует снизить к допустимому уровню. 

Пусть Гах; — максимально допустимая скорость частицы в измерении /. Тогда 
скорость частицы предлагается менять таким образом: 


у, (1+1), если ‚(1+0 =7, 
у, (1+1) = р | 


тах, / 2 


ах,/ > 
если у, (1+0 >И, 


ах, / > 


где значение У,, рассчитывается в соответствии с Без РЗО или фезё РЗО. 
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Величина Га», является довольно важной, поскольку она контролирует общие 
масштабы поиска. Большие значения заставляют мультиагентный метод исследовать 
пространство задачи более тщательно. Если значение Г»; является небольшим, то 
возможны ситуации зацикливания в локальных участках поиска, кроме того, при не- 
больших значениях Г»; увеличивается количество итераций, необходимых для дости- 
жения оптимума. 

Исходя из вышесказанного очевидно, что в разработанном методе величину Гнах; 
предлагается подбирать таким образом, чтобы балансировать между быстрым и медлен- 
ным передвижениями частиц и между фиксацией поиска и общим исследованием про- 
странства. Предлагается рассчитывать Га», таким образом: 

ее т) — 6(х 


тах,/ — Жо» 

ГДЕ Хтах; И Хти;— величины, которые, соответственно, снизу и сверху ограничивают 
диапазон значений измерения /, адеЕе (0,1]. Очевидно, что величина д является проб- 
лемно-зависимой. 

Отметим два важных свойства мультиагентного метода оптимизации с адаптив- 
ными параметрами: 

1) величина Глх; ограничивает не пространство поиска, внутри которого пере- 
двигаются частицы, а только скорости частиц (точнее, диапазон их изменения за од- 
ну итерацию); 

2) максимальная скорость частиц определяется для каждого измерения отдельно 
и обусловливается только его физической природой. 

Исходя из предложенного подхода можно отметить, что ограничение скорости 
частицы может вызвать изменение направления ее передвижения, что может как по- 
зитивно, так и негативно отражаться на поиске оптимума. Этот эффект изображен на 
рис. 1, где х{# + 1) — это позиция 1-й частицы без ограничений на скорость. А величина 


х ‚(1+1) уже отображает позицию частицы после того, как скорость стала равной Га»... 


Нежелательный эффект может возникнуть в случае, если скорости всех частиц ста- 
нут равными Га». Для его предотвращения можно предложить с количеством ите- 
раций уменьшать значение Глх;. Тогда предложенный мультиагентный метод будет 
последовательно сужать масштабы исследования пространства поиска, что делает оп- 
тимизационный процесс более эффективным. 

Тогда мультиагентный метод оптимизации с адаптивными параметрами может 
быть описан в виде последовательности следующих этапов. 


№2 


хкН-1) 
УХа+Т) 


\ 
А 


У (НТ) 


Рисунок 1 — Эффект использования ограничения скорости частицы 
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Этап 1. Инициализация. На этапе инициализации создается множество агентов 
А (14| = п, ), определяющих точки в И. -мерном пространстве поиска. 


Также на этапе инициализации в соответствии с предложенным подходом необ- 
ходимо рассчитать максимально допустимые скорости частиц: 


Г; =0(х —х 


тах, / 
Этап 2. Вычисление значений целевой функции и определение наилучших позиций. 
Данный этап может быть выполнен в последовательности, приведенной выше. 
Этап 3. Обновление скоростей и позиций частиц. Для обновления скорости час- 
тицы предлагается использовать формулу: 


у, (#+1), если у, (1+1) =7, 
у, ((+0=1° Е 
у 


тах, / 2 


тах, / шш, / ) ы 


ах, / > 


если У, ((+0>Г, 


ах, / * 

Для выполнения этапа обновления скоростей и позиций частиц в предложенном 
мультиагентном методе предлагается использовать приведенный выше алгоритм вы- 
полнения аналогичного этапа в РЗО-методе. 

Кроме того на данном этапе необходимо уменьшить величину 6 д= д”, где де (0; 1). 
Уменьшение величины д приводит к тому, что на следующей итерации будет умень- 
шена и величина Га». 

Этап 4. Проверка критериев окончания поиска. В случае, если критерии оконча- 
ния поиска удовлетворены, выполняется останов работы метода, в противном случае — 
переход к этапу вычисления значений целевой функции. 

Разработанный мультиагентный метод оптимизации с адаптивными параметра- 
ми путем управления изменением скорости частиц исключает нежелательный эффект 
слишком быстрого возрастания скоростей частиц, что позволяет улучшить оптимиза- 
ционный процесс за счет уменьшения времени поиска оптимального решения. 


Эксперименты и результаты 


Предложенный мультиагентный метод оптимизации с адаптивными параметрами 
был программно реализован в среде пакета МаЧаь 7.0. 

Используя созданное программное обеспечение, проводились эксперименты по 
поиску оптимального значения функций с целью исследования эффективности предло- 
женного метода по сравнению с известным Р$ЗО-методом. Результаты экспериментов 
с использованием тестовых функций приведены в табл. 1. 


Таблица 1 — Тестовые функции и найденные для них оптимумы 


Оптимум, полученный 
с помощью 
предложенного метода 


Оптимум, полученный 


Функция Фтимум с помощью РЗО-метода 


Ш 


д=У |, 5105) +0А-х, пу: х, = 0,004; х, = 0,0009; 
а | 4 0;0=0. | 200 д.02 
Н5 <х, <5; —5<х.<5. 
12 =20+е зон 0,21 5+ 

‚. > т хи = 0,000299; х = 000005; 
|. оч 5совблт ] а. (0;0) = 0. я = Е 2 —0:000807. 


п=2 
5 <х,<5; —5<х.<5. 
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Главной особенностью предложенных тестовых функций является их полимодаль- 
ность, то есть большое количество локальных оптимумов, что позволяет проверить, 
насколько склонен предложенный метод к зацикливанию в локальных оптимумах. 

Графики тестовых функции представлены на рис. 2 и 3. 


5) 


Рисунок 2 — График тестовой функции Рисунок 3 — График тестовой функции р 


Результаты работы разработанного программного обеспечения для соответствую- 
щих функций изображены на рис. 4 и 5. 


0.8 
——_ Р5О-метод 
06 — — — мультиагентный метод с адаптивными параметрами 
= 
в 0.4 
№ 
0.2 
0 
Номер итерации 
Рисунок 4 — Результаты работы для тестовой функции Л 
5 
Р50О-метод 
4 — — Мультиагентный метод с адаптивными параметрами 


О 10 20 30 40 50 50 70 80 90 100 
Номер итерации 


Рисунок 5 — Результаты работы для тестовой функции р 
Исходя из приведенных результатов можно сделать вывод, что предложенный ме- 
тод позволяет найти более приемлемые решения по сравнению с известным РЗО-мето- 


дом, основанном на мультиагентном подходе. Кроме того, процесс поиска оптимального 
значения происходит быстрее по сравнению с базовым методом. 
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Заключение 


В работе решена актуальная задача ускорения работы мультиагентного метода. 

Научная новизна работы заключается в том, что предложен мультиагентный ме- 
тод оптимизации с адаптивными параметрами, который использует процедуры управ- 
ления изменением скорости перемещения агентов в пространстве поиска, позволяющие 
исключить нежелательный эффект слишком быстрого возрастания скоростей частиц, 
что повышает эффективность оптимизационного процесса за счет уменьшения време- 
ни поиска оптимального решения. 

Практическая ценность результатов работы заключается в том, что разработано 
программное обеспечение, реализующее предложенный мультиагентный метод опти- 
мизации. 
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А.О. Олйник 

Мультиагентний метод оптим1зацй з адаптивними параметрами 

Дослиджено мультиагентний шдх!д до оптим1заци багатовимрних нелийних функщй. Розроблено 
мультиагентний метод оптим1заци з адаптивними параметрами. Проведено експерименти по пошуку 
глобального оптимуму багатовим!рно! функцй на основ! запропонованого методу. 


А.А. ОШпуЕ 

Тве Мшбасепе Орйпи7айоп Метод уйв Адарйуе Рагатеет$ 

'ТБе арргоасН ю орбпиханоп оЁ пи 4итеп$1опа] попПпеаг Влпсйоп$ 15 шуезнзаеа. Тре пивасет( орйпихайоп 
те#фо4 у адарбуе рагатеегз 15 Чеуеюреа. Ехрегитеп оп зеагсН оЁ >1офа] орйтит оЁ пиАнтепзюпа! 
КлпсНоп Базе4 оп Ше оЁеге4 шео4 аге таде. 


Статья поступила в редакцию 28.09.2010. 


90 «Искусственный интеллект» 172011 


